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We study the Shannon entropy of the probability distribution resulting from the ground-state 
wave function of a one-dimensional quantum model. This entropy is related to the entanglement 
entropy of a Rokhsar-Kivelson-type wave function built from the corresponding two-dimensional 
classical model. In both critical and massive cases, we observe that it is composed of an extensive 
part proportional to the length of the system and a subleading universal constant So. In c = 1 
critical systems (Tomonaga-Luttinger liquids), we find that So is a simple function of the boson 
compactification radius. This finding is based on a field-theoretical analysis of the Dyson-Gaudin 
gas related to dimer and Calogero-Sutherland models. We also performed numerical demonstrations 
in the dimer models and the spin-1/2 XXZ chain. In a massive (crystal) phase, So is related to the 
ground-state degeneracy. We also examine this entropy in the Ising chain in a transverse field as an 
example showing a c = 1/2 critical point. 

PACS numbers: 05.50. +q, 75.10.Pq, 03.67.Mn 
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I. INTRODUCTION 

There has been growing interest in quantifying entan- 
glement in extended quantum systems to detect noru 
trivial correlations existing in many-body ground states.EJ 
A useful measure of entanglement is the von Neumann 
entanglement entropy S ,VN (A) := — Tr pA^ogpA, defined 
from the reduced density matrix pa of a subsystem A. Its 
novelty lies in its universal behavior reflecting the long- 
distance nature of the system. In one-dimensional criti- 
cal systems, for instance, the entanglement entropj 
long interval of length I shows a universal scaling 
S' VN (A) ~ | log£ + const., where c is the central charge 
of the conformal field theory (CFT) describing the long 
distance correlations. Possible further information of 
CFT caapbje encoded in a multi-interval entanglement 
entropy.Qul3 As another example, the existence of topo- 
logical order in gapped systems can be detected by 
measuring a constant contribution to the entanglement 
entropyall3 (with, recent applications to^^tioiiaLfluan- 
tum Hall stately and Z 2 spin liquidsHjIlaOO) . 

Here we introduce an apparently different entropy 
as follows. Consider a one-dimensional (ID) quantum 
model and its ground state \g). If one chooses an or- 
thogonal basis {(£)} of the Hilbert space, one gets a set 



of probabilities pi := 
tropy can be defined: 



', from which a Shannon en- 



s' := - ^Pi log Pi 



(1) 



Note that this entropy depends on the choice of basis. 
For a f7(l)-symmetric model with the conservation of the 
particle number or the magnetization, we use the local 
particle occupations {nj} or magnetizations {c|} to de- 
fine the basis. The entropy S is small when the wave 
function \g) is dominated by a particular crystal state 



| z o ) - It becomes larger as more basis states contribute 
to the wave function due to quantum fluctuations. Thus, 
this entropy quantifies quantum fluctuations or entangle- 
ment occurring in the given basis. Like other entangle- 
ment measures, we will see that the scaling of this entropy 
is controlled by the essential long-distance nature of the 
system. We note that a similar enteopry also appears in 
the context of dynamical systems jLM3 where it is used 
to quantify chaos.E2 

At the same time, this entropy has an interpretation 
as the entanglement entropy of a two-dimensional (2D) 
quantum state, as we describe in detail in Sec. 0. The ba- 
sic idea goes as follows. A ID quantum Hamiltonian on 
a ring is related to a 2D classical model on a cylinder in 
the transfer matrix formalism- Then a Rokhsar-Kivelson 
(RK) -type wave functionESO |RK) can be constructed 
from this 2D classical model. One can show that in the 
limit of a long cylinder, the entropy S defined in Eq. (Q) is 
precisely the entanglement entropy S VN (A) of the 2D RK 
state |RK) for a half cylinder A shown in Fig. Q. More 
specifically, each probability pi becomes an eigenvalue 
of the reduced density matrix pa- This correspondence 
allows us to apply a simpler ID picture to study the uni- 
versal behavior of entanglement entropy in 2D pitkjal 
systems, a very active subject in recent literature E3EMfl 

The purpose of the present paper is to unveil the 
generic scaling properties of the Shannon entropy (|l|) of 
ID ground states \g), as a function of the ring length L. 
This amounts to studying the entanglement entropy S VN 
of 2D RK states |RK) defined on a cylinder of circumfer- 
ence L. In both critical and massive systems, we observe 
that for large L, the entropy is composed of an extensive 
part proportional to L and a subleading constant: 



S{L) = aL + S Q + o(l). 



(2) 



The extensive part aL simply reflects the fact that a 
generic wave function \g) spreads over an exponentially 



2 



large number of microscopic configurations. In terms of 
the 2D entanglement entropy, this can be interpreted as 
a boundary contribution. Here, we are however inter- 
ested in the subleading constant So- As we will see, So 
is universal and is determined by the basic properties of 
critical or massive systems. 

Our primary interest lies in the situation where a ID 
quantum or 2D classical system (used to build a RK 
state) is described by a c = 1 massless bosonic field 
theory (Tomonaga-Luttinger liquid; TLL) with the Eu- 
clidean action: 



1 

8^ 



dxdy [{d^f + (d y <pf 



(3) 



Here the bosonic field is compactified on a circle: 4> = 
</> + 2irR. The boson compactification radius R is an im- 
portant scale- invariant number which controls-the power- 
law behavior of various physical quantities £a We find 
that So is given by a simple function of the radius: 



So 



log*--. 



(4) 



We present detailed analyses to establish this result. 
In Sea S H I we study the Dyson-Gaudin Coulomb gas 
model,Effi3 which gives probabilities {pi} for the dimer 
model on the hexagpnal lattice and the Calogero- 



III C, we 



Sutherland modelEZEa In particular, in Sec. 
analytically derive Eq. (|j) using a free field representa- 
tion of the gas model. In Sec. IV, we numerically demon- 
strate the same result in the spin-1/2 XXZ chain. 

At a certain value of R, the system undergoes a phase 
tran sition to a massive crystal phase. As we show in 
Sec. [Ill D| and Sec. |v|, in the massive phase breaking a 
symmetry, So is related to the ground-state degeneracy d: 



S = log d. 



(5) 



At the transition point, we observe a jump in So, though 
it is slightly obscured due to finite-size effects. 

As an example showing the c = 1/2 criticality, in 
Sec. 0, we study an Ising chain in a transverse field (63). 
We calculate the entropies in a z and a x bases, corre- 
sponding to the RK states built from an eight-vertex 
model and a 2D Ising model, respectively. The extracted 
constant S^ ~ -0.4387 and S { x) ~ 0.2544 at the crit- 
ical point might be generic constants characterizing the 
c = 1/2 CFT, although we do not have any analytical 
derivation of these numbers. In the symmetry-broken 

(z) 

phase, the constant takes Sq = — log 2 in the a z basis, 
which is interpreted as a manifestation of Z2 topological 
order in the eight-vertex RK state. 

As a related quantity, in Sec. VI, we study the scaling 
of the largest probability po '■= m&xpi. This maximum is 
attained by crystal states |zo), e.g., by Neel states |flTI 
. . .) and lltlt • ■ ■) f° r the XXZ chain in zero magnetic 
field. Very similarly to the entropy S(L), the logarithm 
of po has a dominant linear contribution followed by a 
subleading constant: 



For c = 1 critical systems, our numerical results in the 
Dyson-Gaudin gas and the XXZ chain show that 



7 = logi?. 



(7) 



This result, together with Eq. gives a simple and 
universal way to determine the radius R from a ground 
state wave function. 

Here we comment on closely related works. Using 
boundary CFT, Hsu et alx3 have also studied the en- 
tanglement entropy of 2D RK wave functions for a half 
cylinder. Their prediction for the constant So differs from 
ours, although it matches our calculation for a different 
constant 7 appearing in — logpo- A quantity similp; to 
Po has also been studied by Campos Venuti et a/.c3 in 
the context of fidelity, in agreement with our result (Q) . 



II. CORRESPONDENCE BETWEEN SHANNON 
AND ENTANGLEMENT ENTROPIES 

In this section, we formulate the connection between 
the entanglement entropy of Rokhsar-Kivelson states and 
the Shannon entropy of the ground states of ID quantum 
models. 



A. Generalized Rokhsar-Kivelson states 

We start from a discrete classical (spin) model on a 
2D lattice defined by Boltzmann weights e~ E ^ for mi- 
croscopic configurations c of the system. The partition 
function is given by 



Z 



E< 



-E(c 



(8) 



The Hilbert space of a 2D quantum system is constructed 
by associating a basis state |c) to each classical config- 
uration c. Then, one can define a generalized Rokhsar- 
Kivelson (RK) wave function as the linear combination 
of all the basis states |c) with amplitudes given-nte, the 
square roots of the classical Boltzmann weights:! 



|RK) 



1 



E 



,-h E ( c )\ 



(9) 



logpo = olL + 7 + o(l). 



(6) 



The RK state shares the same correlations with the orig- 
inal classical model, as long as one focuses on the diag- 
onal correlations in the {|c)}-basis. This type of state 
was first introduced as the ground state of the quantum 
dimer model (QDM) on the square lattice, where {c} 
are fully packed dimer coverings of the lattice. Later, 
the same type of states have been studied_£or different 
lattices (hexagonal,E3 triangular,!^, kagomeEj, etc.) and 
for a modified dimer modelHa The RK states for the 
eight-vextex model on the square lattice have also been 
studiedEil 
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FIG. 1: Two dimensional system with a cylinder geometry, 
divided into upper and lower parts, A and B. 



B. Schmidt decomposition and entanglement 
entropy 
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FIG. 2: (color online) Spatial division of the square latti ce into 



regions A and B, corresponding to the procedure in Sec. II B 1 . 
The system is defined on a cylinder, and is periodic in the 
horizontal direction. 



We divide the system into two parts, A and B, as in 
Fig. U The reduced density matrix pa is obtained from 
a state |RK) by tracing out the degrees of freedom in B: 



p A = Tr |RK)(RK|. 



(10) 



We are interested in the (Von Neumann) entanglement 
entropy of A: 



S VN (A) = -Tr p A log pa- 



(11) 



Here we show that the calculation of S" VN (A) for a RK 
state @ can be recast as a fully classical calculation, 
provided that the boundary between A and B satisfies 
certain geometrical conditions. This is done by deriving 
a Schmidt decomposition of the RK state @ . 



1. Case with local constraints 

We first consider the case where the classical model 
contains certain local constraints. For simplicity, we as- 
sume that the system consists of Ising variables Oj sit- 
ting on the bonds of the square lattice as in Fig. | (The 
same argument also applies to a model on the hexago- 
nal lattice.). Around each site, the four Ising variables 
Uj :'s satisfy the following constraint: if three of them are 
specified, the last one is uniquely determined. Dimer, 
six-vertex, and loop models satisfy this condition. For 
example, in the case of a dimer model, we assign aj = ±1 
to the presence/absence of a dimer on the bond j. Then, 
there is strictly one bond with a 3 ■ = +1 emanating from 
each site. 

We define the system on a cylinder, and divide it into A 
and B, as shown in Fig. ||. Here, all spins at the boundary 
belong to A. The spin configurations a and b inside A and 
B must agree with the configuration i at the boundary. 
Let Ef- (£P) be a set of such a's (6's). Thanks to the 
local constraints, such sets share no common element: 



£? n£<? = (CI = A, B; i^i'). 



(12) 



We assume that the classical model E(c) contains only 
local interactions involving four bonds around each site. 
The energy can then be decomposed into two parts: 



E(c) =E A {a,i) + E B (b,i). 



(13) 



The first term corresponds to all the interactions among 
spins in A. The second corresponds to interactions 
among spins in B and the boundary region. The im- 
portant point is the absence of direct interaction between 
spins inside A and B. Thanks to this property, the Boltz- 
mann weight of the configuration c = (a,i,b) factorizes 
into two parts, which allows us to rewrite Eq. (||) as fol- 
lows: 



e-i EA ^\a,i] 
J2 e~^ E ^\b) 



be£f 



(14) 



We define normalized RK states (boundary-dependent) 
on A and B as 



irk; 



IRK, 



_1_ £ e-i^-Oio, 



a££f 



_J_ £ e -i*»(M)| 6> 



be£f 



with 



(15a) 
(15b) 

E^,i) (n = 4,B).(15c) 



Then we arrive at the Schmidt decomposition 

|RK) =Y,^ |RKf >|RKf >, with Pl := ?L0L. 



(16) 
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FIG. 3: (color online) Spatial division of the square lattice, 

Each spin of 



1IB2 



corresponding to the procedure in Sec. 
the original model is replaced by a pair of spins. An infinitely 
strong ferromagnetic interaction ensures that the two spins 
are always in the same state. The (new) lattice is divided in 
two regions, A and B, and the original sites which overlap 
with both regions are called "boundary" sites. The RK state 
admits a simple Schmidt decomposition (see text) if the only 
interactions between A and B take place inside the boundary 
region. For general short ranged interactions (not necessarily 
first neighbor), this condition can be achieved by choosing a 
sufficiently "fattened" boundary region. 



Here, the mutual orthogonality (RKf|RKp) = 6u> is 
guaranteed by Eq. (|l^). The reduced density matrix pn 
(with CI = A or B) is thenEl 



= 5>|RKf)(RI< 



(17) 



freedom lives on every site. We again assume that the 
energy E(c) consists only of interactions between nearest- 
neighbor spins. As illustrated in Fig.||, each spin is dupli- 
cated and an infinitely strong "ferromagnetic" interaction 
is added so that the two copies of the original spin are 
always in the same state (no spurious degrees of freedom 
are introduced). Then the two regions A and B are in- 
troduced in such a way that all the spin-spin interactions 
in E(c) take place inside A or B. In other words, the 
only allowed couplings between A and B are the infinitely 
strong "ferromagnetic" interactions between copies of the 
same physical spin. If one prefers to think in terms of 
the original spins only (not duplicated), this amounts to 
saying that the regions A and B are overlapping around 
their boundary. 

In this setup, each state |c) can be labeled in the fol- 
lowing way: 



\a, i) ® \b, i). 



(19) 



Here, each original spin lying at the boundary is effec- 
tively "split" and has one copy in |o,i) and the other 
in \b,i). The Schmidt decomposition ( |l6| ) is then con- 
structed using the following states: 



|RKf) 
with 



1 _i 
-== > e 2 



U3, l) 



(20a) 



E< 



-E n (ui,i) 



(SI = A, B). (20b) 



Here, the difference from Eq. ( |l5j ) is the presence of i 
inside the ket \w, i) for both SI = A and B, which ensures 
the mutual orthogonality of {|RKf )}. 



Therefore, we get 

5 VN (i) = -^ K logft. 



(18) 



The entanglement entropy can thus be computed from 
the classical partition functions Z^Zf with boundary 
spins fixed in a state i. Similar formulations were used in 
Refs. [T^ , ^4] , p!5|Jl^ for exact /perturbative/numerical cal- 
culations of entanglement entropy in toric code, quantum 
eight-vertex, and quantum dimer models. 



2. General case 

The above discussion relies on the presence of local 
constraints. Without them, any configuration b is al- 
lowed i n B irrespective of i, and states {|RKf )} defined 
in Eq. (15b) are not mutually orthogonal in general. Even 
in such a general case, one can slightly modify the model 
so that the entanglement entropy can be computed in the 
same formulation. 

For simplicity, we assume that the classical model is 
defined on the square lattice, and a spin-S' degrees of 



C. Transfer matrix calculation of the reduced 
density matrix spectrum 

In the previous subsection, the spectrum {pi} of the 
reduced density matrix has been expressed in terms of 
the classical partition functions with spins fixed in the 
boundary region. In the cylindrical geometry of Fig. |l| 
with circumference L and height 2h, we can relate this 
spectrum to the ground state of a ID quantum spin model 
using the transfer matrix formalism. 

Corresponding to the 2D classical model, we introduce 
the transfer matrix T in the upward direction in such 
a way that it connects spin configurations on neighbor- 
ing "rings" winding around the cylinder (e.g., a\ and 
shown in Figs. || and ||). The classical partition function 
(||) is then expressed as 

z = E E E 

a h -i,...,ai i bi,...,bh,—i 

(o h |T|o h _i) . . .(oalTloiXoilTI*) (22) 
x(»|r|6i)(6i|T|62>...(6 h _i|r|6 h ) (23) 
= (a h \T 2h \b h ). (24) 



5 



Here, the spin configurations, ah and bh, at the top and 
bottom edges of the cylinder are fixed. In this setup, the 
classical probability to find a given configuration i on the 
ring (boundary) is 



1 



Pi 



\a h \T h \i)(t\T h \b h ). 



(25) 



S(J3 = 2) = S VN (A) to arbitrary (3 > 0. It can be com- 
puted through the standard thermodynamical relation 
S(/3)^ (l-0d )Z(f3). This formulation will be useful in 



Sec. [II. We note that a similar extension of entanglement 



entropy has also been discussed by Li and Haldane.1 



We now consider the limit of a long cylinder h L 
so that only the dominant eigenvector \g) of T (with 
the largest eigenvalue too) contributes. Using T h ~ 
™ft\g){g\, we get: 



Pi 



\(i\g)f 



(26) 



If the transfer matrix is related to ID quantum Hamilto- 
nian H via T ~ e~ T7i (with t being a small time interval), 



\g) is the ground state of H. Then, Eq. fl2q) means that 
the complete spectrum of the reduced density matrix pa 
is given by the elements of the ground state vector \g). 
In spite of its simplicity, to our knowledge, this result has 
not been reported previously in the literature. 

In the rest of the paper, we will make an extensive use 
of this property to calculate the entanglement entropy 
(|l~8|). We study several RK states (dimer, vertex, and 
Ising models) defined on the infinite cylinder for relatively 
large values of L (maximally, L = 48 for the dimer models 
and L = 32 in the six- vertex models with no field). 



D. Thermodynamic extension 

The spectrum {pi} of the reduced density matrix pa for 
a RK state has a simple interpretation in terms of bound- 
ary free energy of the classical model. Using Eq. (|l6|), we 
get 



F B - F 



(27) 



where Ff- :— — \ogZf- (resp. F^) is the free energy of the 
subsystem A (resp. B) with its boundary with B (resp. 
A) fixed in the configuration i. Also, F = — \ogZ is the 
free energy of the whole system, without any constraint 
on the spin configuration at the boundary between A and 
B. We now identify the r.h.s. of Eq. ^2^ ) as an effective 
energy 2E(i) for the boundary spins OThen the entan- 
glement entropy 5 VN (.A) in Eq. ( |i~8| ) can be interpreted as 
the "thermal" entropy for the boundary spins. We push 
further this thermodynamic interpretation of S (A) by 
introducing a parameter in order to modify the proba- 
bilities pi: 



Pi 



-2E(i) 



Pi 08) 



1 



Z(J3) 



-0E(i) 



(28) 



where Z(j3) = ^ e - ' 3 ^ 4 ) is a normalization factor (with 
Z(f3 = 2) = 1). Here, (3 plays the role of an effec- 
tive inverse temperature for the boundary spins (but a 
priori not for the bulk of the classical system defined 
by Eq. (§h). This allows us to generalize the entropy 



III. FROM CRITICAL DIMER RK STATES TO 
DYSON-GAUDIN COULOMB GAS 

In this section, we study critical dimer RK states on 
bipartite lattices and a related ID classical gas model. 
Using the formulation of Sec. ||, we compute the en- 
tanglement entropies of the dimer RK states using the 
dominant eigenvectors of the transfer matrices. As de- 
scribed in Appendices and [C], the transfer matrices of 
the dimer models can be expressed as free fermion Hamil- 
tonians, and their dominant eigenvectors are Slater de- 
terminants. In particular, for the hexagonal lattice dimer 
model, the resulting probabilities {pi} coincides with the 
Boltzmann weights of a ID lattice gas interacting via 
a repulsive logarithmic potential (Dyson-Gaudin gas) at 
the inverse temperature £L.= 2. This gas model is also 
related to the discretizea£3 ground state wave function 
of the Calogero-Sutherland model. Therefore, the entan- 
glement entropy of the dimer RK state and the Shannon 
entropy of the discretized Calogero-Sutherland ground- 
state wave function coincide with the thermal e ntrop y of 
the g as model. As we will demonstrate in Sees. [II A] and 
HI B , this entropy contains a no n-extensive constant con- 
tribution, Eq. ([I). In Sec. Ill C , we derive it analytically 
using the free field representation of the gas model in the 
continuum limit. 



A. Critical dimer RK states 

We start from the RK states (^) constructed from the 
dimer models on the hexagonal and square lattices. The 
energy E(c) takes zero for any fully packed dimer con- 
figuration c and infinity otherwise. These RK states are 
the ground states of quantum dimat,p,odels at special 
points and are know to be critical.caEJ By associating 
an Ising variable Oj = ±1 with the dimer occupation on 
each bond j, w e can adopt the formulation presented in 
Sees. [I B 1 and II C to compute the half-cylinder entan- 
glement entropy S' vN . 

As described in Appendices § and ^, the transfer ma- 
trices T of the dimer models are expressed as free fermion 
Hamiltonians. Their dominant eigenvectors \g) are Slater 
determinants (Fermi sea), and the weight pi = \{i\g)^ 
can be computed by evaluating a determinant. Let us 
focus on the hexagonal case, where a simple expression 
for {pi} is available. We here assume that L is a multiple 
of 6 for the sake of simplicity. A generic configuration i 
of the boundary will be given by a number n G {0, . . . , L} 
of fermions living on the vertical edges of the boundary 



and their positions: 



< aci < . . . < a n < L - 1. 



(29) 



It is also shown in Appendix [Bjthat, in the limit /i> L, 
the only nonzero probabilities correspond to n = 2L/3 
fermions and are given by a Vandermonde determinant, 
which simplifies into: 



^ n 



4 sin 



i<j<i'<« 



-(ptj - a y 



(30) 



This equation is invariant under n — > L — n. Therefore 
it is easier to compute it with n = L/3 fermions, instead 
of n = 2L/3. As we will see in the next section, these 
probabilities coincide with the Boltzmann weights of a 
ID lattice gas. Note that the present calculation is done 
without any constraint on the configurations, ah and bh, 
at the top and bottom of the cylinder. The number n 
of fermions leading to nonzero probabilities can be con- 
trolled by imposing certain boundary conditions. 

Figure [| shows the scaling of the entanglement entropy 
S(L) for both the hexagonal and square cases. Here, in 
the square case, the weight pi is co mpute d by numerically 
evaluating the determinant in Eq. ( C1S| ). In the hexago- 
nal case, we examine different fermion densities p = n/L 
(or equivalently flux/ winding sectors in the dimer lan- 
guage). In every case, the scaling of S(L) appears to be 
approximately linear in L. The slope is non-universal, as 
expected, and depends on the details of the system. The 
most interesting result is the existence of a finite constant 
So- By fitting the large-L values by S(L) = aL+So+b/L, 
we find So = —0.500 ± 0.002 in all cases. These results 
suggest that Sq = —1/2 is a universal number for this 
family of RK wave functions. 



B. Dyson- Gaudin gas on a circle 

We consider a system of n charges Q = — 1 , living on a 
periodic one-dimensional lattice with L sites as in Fig. [|. 
These charges interact via a 2D Coulomb repulsive po- 
tential equal to minus the logarithm of their mutual dis- 
tance. For convenience, we add a uniform background of 
L charges Q = +1/2 located on each site of the lattice. 
A configuration of the system will be determined by the 
positions 1 < ai < a-2 < . . . < a n < L of the charges. 
Two charges cannot occupy the same site. The energy of 
a given configuration i is: 



- E 

l<J<j'<" 



log 



e £ 



e i 



logi, (31) 



where the second term comes from the interaction of each 
charge with the background. ThisJiiay be viewed as a dis- 
cretized vepion of the Dyson gas,t3 and has been studied 
by Gaudi nEa (h ence we name it the Dyson-Gaudin gas). 
As in Sec. Ill A , the "particle-hole" symmetry n — > L — n 
is worth noticing. The classical ID partition function is 



Dimers on the hexagonal and square lattices 
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FIG. 4: (color online) Entanglement entropy of RK states 
corresponding to dimer models on the hexagonal lattice (with 
fermion density p = 1/4, 1/3, 1/2) and the square lattice 
(with p — 1/2). In all cases, the entropy scales as S ~ aL + 
So + b/L with So = -0.500(2). The inset shows So, the 
subleading constant computed this time from a linear fit S ~ 
a L + S'o on the interval [12, L] as a function of 1/L. The 
convergence towards —0.500(1) can be seen. 




FIG. 5: System of n = 4 charges Q = — 1 on a circle with 
L = 12 sites. The left shows the configuration with highest 
probability (degeneracy: 3) and the right shows a random 
configuration. On each site, a "background" charge Q = +1/2 
is added. 



given by Z{8) = J2 { < 



-PE{i) 



for an inverse temperature 



0, and the classical entropy can be calculated from it: 



S((3) 



-4 



]og Z(J3). 



(32) 



Gaudmll has evaluated the partition function exactly in 
the special cases where 3 — 2A with A an integer and 
A < L/(n- 1): 



ZW(2A) = 



(nA)! 



n!jL n(A-l)( A l)r 



(33) 



The probabilities {pi} in Eq. ( |3C| ) calculated for the 
hexagonal lattice dimer model coincides exactly with the 
Boltzmann weights of the Dyson-Gaudin gas model with 

13 = 2: Pl 
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Hence the entanglement entropy for this dimer model 
corresponds to the thermal entropy S(ft = 2) of the gas 
model. Note that Z(ft = 2) = 1 if we set A = 1 in 
Eq. (33) because pi i n Eq , (pQ ) are already normalized. 



Dimers on the hexagonal lattice : sub-leading constant S in the entropy 



In the spirit of Sec. II D , we can generalize the entan- 
glement entropy by changing the inverse temperature ft. 
We define Pi {ft) = pf /2 /Z(ft) = e~^ Ei /Z{ft), and the 
associated Shannon entropy 



(34) 



coincides with the thermal entro py (|32| ) of the gas model. 
Notice that Gaudin's solution (|33| ) cannot be used to 
compute the entropy S((3) because it is valid only for 
special values of ft. Instead, we compute S((3) numeri- 
cally by explicitly summing over all the configurations. 

Here we mention the connections with other models. 
The Dyson gas model emerges in the weights |(i|g)| 2 
of the Jastrow-type ground-stata-wajre function of the 
Calogero-Sutherland (CS) model£3u3 The inverse tem- 
perature ft of the gas model is related to the coupling 
constant of the CS model. The CS model is described as 
a Tomonaga-Luttinger liquid at low energies. According 
to a spectral analysis of the CS modelJ^J the boson radius 
R is related to ft as 

i?= vf- (35) 

This allows us to control R simply by changing ft. This 
relation will also be justified from a different vie*w>«int 



in the next section. The Haldane-Shastry modelE§E3 is 
a discretized version of the CS model, and its Jastrow 
ground state (Gutzwiller-projected Fermi wave function) 
has the weights exactly obeying the Dyson-Gaudin gas 
with ft = 4. This model has R = s/2 because of the 
SU(2) symmetry, in consistency with Eq. (|35|). Note also 
that the same wave function is known to be an extremely 
good ansatz for the ground state of the Heisenberg chain. 

Now we analyze the thermal entropy S(ft) of the gas 
model. We extract the non-extensive constant contribu- 
tion So as in Fig. [|, and plot it as a function of R = y/ ft/2 
in Fig. ||. We find that the data agrees well with a simple 
relation 



So =logi? 



1 



(36) 



This expression is derived analytically in the next section. 
It should be noted that the subleading constant (|36| ) is 
increasing with ft, contrary to the total entropy, which 
is decreasing with ft, as should be the case in classical 
thermodynamics. 



C. Free bosonic field 

In this section, we obtain the expression of the entropy 
So using a field theoretical approach. Our goal is to ob- 
tain a continuous expression for the partition function of 



Subleading constant (up to L=39) 
log[R]-1/2 



0.6 0.8 



1.2 1.4 



FIG. 6: (color online) The subleading constant So in the ther- 
mal entropy of the Dyson-Gaudin gas with density p = 1/3, 
as a function of R = \J ft/2. At ft = 2, this coincides with 
the constant part in the entanglement entropy of the hexag- 
onal dimer RK state shown in Fig. ^. The data well obeys 
Eq. fjjj). 



the gas studied in the last section and deduce from it the 
expression of the entropy. 

We consider a continuous distribution, p{9), of electric 
charges on the unit circle. The expression of the electro- 
static energy Eq. (|3l|) is given by: 



E 



p{0)d6 / p(9')d6' log 



2 sin 



■(37) 



We define a field (f)(9) measuring the amount of charge in 
the interval [0,9] in units of 27r: 



(j}{9) = 2ir / p(a)da 
Jo 



(38) 



By performing partial integrations twice, the energy is 
rewritten as 



(39) 



Since the functional integration is over p = <f>' , the 
zero mode of (j> is unphysical and should be discarded 
(it can be removed by adding an appropriate constant to 
Eq. @)). By expanding the field cj) over modes 




o_ \ ^ im9 , — —irn8\ 

2n [x m e + x m e ) 



the energy Eq. (|39| ) reduces to 



E 



E 



m\x m \ 



(40) 



(41) 



Integrating the Boltzmann weight e /3E over the modes 
dx m dx m , we obtain the partition function Z sp h ro of the 
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'as: 



^sphere 



n 



2n/3m 



(42) 



To find possible universal contributions to the corre- 
sponding free energy, this expression of course needs to be 
regularized. Following Nahm,c3 we regularize the mea- 
sure dx m dx m to take into account the finiteness of the 
number of states. We set |a> m | = p m and take the mea- 
sure to be d(2n[pj n }^ m ^ A ^) where f(x) = 1 in the interval 
[0, 1] and decreases to /(oo) = sufficiently fast. We ob- 
tain : 



^sphere 



n 



i 



/(m/A) 



T(l + /(m/A)) (43) 



The Euler-MacLaurin formula yields : 



rn—l v 



r(i+/(=J) 



/^ e A/ oo ^(x)dx (1+o(1)) 



gp{x) = lo g r(l + f{x)) - log(7r/3)/(x) (44) 

Here we decompose / into the sum of a step function on 
[0, 1] and a function vanishing on [0, 1] : 



n - 

\ 111 



fi 1 



(45) 

where we have used the Stirling formula on the term cor- 
responding to the step function, and Euler-MacLaurin's 
on the other one. Putting everything together : 



= y/ZeA/o 00 ^lo g r(l+/(x))-log(7r^Ax)/( : c)^ 1 + Q ^ 

(46) 

The regularized partition function is obtained after re- 
moving the exponential factor, A~ aA p A , which can be 
thought as the "extensive" part 



^sphere 



^sphere 



(47) 



This result is simply equal to the C regularization of 
Eq. (|4^). In this derivation, the normalization of the 
integration measure over the modes is adjusted so that 
•^sphere = 1 at (3 — 2, to agree with our microscopic defi- 
nition of the partition function, Z sp h crc (ff) = Yl<iPi > m 
the discrete model. From this expression, we deduce the 
thermal entropy: 



So 



1 



log 



(48) 



in agreement with our numerical result in Fig. |6| 

Alternatively, we can normalize the field 4> differently 
so as to include the inverse temperature f3 in its defini- 
tion. We introduce a radius R and set j3 = 2R 2 . The field 



4> is now defined modulo 2itR. We can extend its range 
of definition by requiring it to be a harmonic function 
on the unit disk fi, and express the energy as a Dirichlet 
integral: 



2 4?r 



dzdz d z 



(49) 



where z = x + iy. This coincides with the action (||) 
except that the range of integration is limited to the unit 
disk Q. | — | 
One can view Eq. ( |49| ) as the action of a closed stringpil 
propagating on a circle of radius R with a Regge slope 
a' = 2. Now, if we include the center of mass (zero mode) 
into the definition of <j>, we see that the regularized mea- 
sure [d<p] becomes invariant under rescaling of the field, 
so that we can also take the field defining the measure to 
be defined modulo 2-kR. The partition function Z sp hcrc of 
the electrostatic gas is obtained by sewing together two 
disks to form a sphere and is given by the partition func- 
tion of a closed string propagating on a circle of radius R. 
Proceeding in this way, the oscillators do not contribute 
and it reduces to the center of mass integral: 



^sphere 



R, 



(50) 



-A r °°lo g (Ax)/(x)dx (1 + o(1)); 0( 



in agreement with Eq. (|47|). 

To understand this result, consider a closed string 
cr, T) propagating in the Euclidean time T. Its parti- 
tion function on the cylinder [0, 27r] x [0, T] with boundary 
fields equal to 4>i t 2 defines the propagator Z{4>\, 4>2)- We 
evaluate the torus partition function Z tom s by taking the 
trace of the propagator over 4> — 4>i — 4>2 ■ 



]Z< 



(51) 



If we decompose the field 4> into the sum of a harmonic 
function equal to 4>i, 4>2 at the two boundaries and a field 
vanishing at the boundaries, we factorize the propagator 
into two pieces: a classical one equal to e"! 6 *)"!*) 
in the limit of large T, and the partition function Zqq 
with Dirichlet boundary conditions. Thus, In the limit 
of large T, Z toms = Z sp h e re^oo- But Z m /Zt oms = 1/R 
is the stationary probability distribution of the center 
of mass of the string diffusing on a circle of radius R 
and the result Eq. ( p0|) follows. This formulation has the 
advantage to be easily generalizable to a closed surface 
of characteristic \'- 



Z, 



(52) 



We expect this formula also to apply to the case of open 
boundary conditions with suitable chiral boundary con- 
ditions. 

Let us mention that Gaudin's partition function (33) 
has the same asymptotic expression as Sphere m Eq. (42) 
if we remove the non-universal extensive part. Indeed, 
by taking n — > oo while keeping p = n/L constant, the 
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Stirling's formula nl 
gives: 



/ 2m (n/e)" applied to Eq. (§|) 



Zi L) (2X) = fi(\) n R, with fx(X) 



3 a-i a a 6 a-i 
A! 



(53) 



Gaudin's formula (|33|) is only valid at integer values of 
A and we cannot use it to evaluate the entropy. How- 
ever, the denominator factor L< A_1 )™ in Eq. (|33| ) may 
be viewed, pas a regularization of the Dyson partition 
function J£3'c3 where L plays the same role as A defined 
above, and the entropy also derives from the Dyson gas 
partition function. 

The compartjfi cation radius of the (continuum) CS 
ground stateE3£3 is known to be R = ^/(A = (3/2) from 
Ref. On the other hand, our numerics on the dis- 
cretized version of the CS wave function indicate that 
the entropy constant is equal (within our numerical ac- 
curacy) to So = log(-\//3/2) — 1/2. From the analytical 
derivations presented above, this entropy constant must 
be related to the boson radius. We therefore conclude 
that the discretized CS state has long distance properties 
described by the same boson compactification radius as 
the original continuum CS wave function. It is interesting 
to notice that this identification has been made through 
the ground-state structure of the CS model, i.e. with- 
out relying pon the spectral properties unlike preceding 
approaches £3 



D. Phase transition toward a crystal state 

The numerical results of Fig. |o| show that for not too 
large values of f3, the non-extensive contribution to the 
entropy is given by So = \ogR — 1/2. As should be clear 
from the previous section, this can only be true if the 
system is described by a massless (but compactified) free 
field. But for sufficiently large f3, the system undergoes a 
transition to a crystal state with spontaneous translation 
symmetry breaking. A simple way to understand that 
such a crystal is expected at large (3 is to notice that 
for (3 — > oo, only the particle configurations i for which 
the (original) probability pi is maximum survive. For a 
particle density 1/d, this selects d periodic configurations 
with equally spaced particles (see the left panel in Fig. ||). 
Adding fluctuations around these regular configurations 
will add extensive contributions to S(L), while keeping 
the subleading constant 



S = log d 



(54) 



stable in a crystal phase with a d-fold spontaneous sym- 
metry breaking. Due to some finite-size effects, it turns 
out that the liquid-crystal transition is easier to see in 
the non-extensive part C of losZ (rather than that of 
S). The data displayed in Fig. |7] are consistent with 



C = 



log R, R<R c = d 
logd, R>R c = d. 



(55) 



Subleading constant in log(Z) : transition towards a cristal state 

Hexagonal : density p=1/3 degeneracy=3 — I — 
Hexagonal : density p=1/2 degeneracy=2 

Square : density p=1/2 degeneracy=1 ■■■^K— ■ 

log(R) 

log(d) 




R=((5/2) 1/ ' ! 

FIG. 7: (color online) Liquid-crystal transitions induced by 
the boundary temperature (3" 1 in the dimer models. Here, the 
subleading constant C in log Z is shown. The critical radius 
is R c = d (/3 C = 2d 2 ), where d is the degeneracy of the ground 
state. The constant C is expected to obey Eq. Jj^). The 
discrepancy slightly after R c is very likely due to finite-size 
effects. 



Using So = (l — <9r) C, we can recover the subleading 
term Sq in the entropy: 



So 



logi?- 1/2 
\ogd 



R < R c 
R > R r 



= d 
d. 



(56) 



It should be noted that the transition is only visible on 
the subleading terms of S and \ogZ. 

The crystallization can also be understood from a free 
field point of view. Let us perturb the action (||) by a 
d-fold symmetry breaking boundary field : 



= / p h d So* cos(dct>)de\ 



I sphi 



(57) 



where the integral is taken over the equator of the sphere. 
In a spin wave approximation, the anomalous dimension 
of the field hd is: 



1 d2 



(58) 



Thus, the perturbation becomes relevant when R > d in 
agreement with our observations (Fig. [?]) . 

This transition is completely analogous to the localiza- 
tion transition of a macroscopic degree of freedom cou- 
pled to a dissipative environment in presence of a peri- 
odic potentials In this context, the inverse temperature 
(3 = 2R 2 is a friction coefficient. A similar transition is 
observed in the XXZ chain (see the following section), 
but at a different value of the compactification radius 
(i? = \/2), compatible with the bulk roughening transi- 
tion of the 6-vertex model. 
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IV. SPIN- 1/2 XXZ CHAIN AND SIX- VERTEX 
RK STATES 

In this section, we consider the Shannon entropy (0) 
defined from the ground state \tp) of the spin-i XXZ 
chain: 



Spin 1/2 XXZ chain : subleading constant in the entropy 



(59) 



This Hamiltonian is related to the transfer matrix, of the 
classical six- vertex model on the square lattice Thus, 
using the argument of Sec. ||, the Shannon entropy S 
here can also be interpreted as the entanglement entropy 
of the RK state built from this vertex model. Since the 
magnetization per site M — \ J2i °f ls a conserved quan- 
tity, we can work in a sector with fixed M. We calculate 
the ground state of Ji for finite periodic chains using 
Lanczos diagonalization (up to L = 32 for M = and 
L = 40 for M = 1/2), and evaluate the Shannon entropy 
S from it. 

We first focus on the c = 1 Tomonaga-Luttinger liquid 
phase extending over a wide region in A > —1. The 
boson radius R depends on A and M (see Refs. £l6] and 
for details). When M = 0, R is related to A via a 
simple relation: 



Exact diag M=0 -\- 
Exact diag M=1/5 
Exact diag M=1/4 ^ 
Exact diag M=1/2 [ 
log[R]-1/2 




0' 



A=-0.6, M-0 
A=0. NM) 
A=0.6. M=0 



R 

FIG. 8: (color online) The subleading constant So in the en- 
tropy S extracted from the critical ground state of the XXZ 
chain (|H^). The examined values of A range from —0.8 to 
1 for M = and from -0.8 to 8 for M = 1/5, 1/4, and 
1/2. The inset shows the fitting of the data with the scaling 
form S = aL + So + b/L. The constant So well obeys the 
proposed universal formula log R — ^ . Close to the isotropic 
point (A = 1 and M = 0) with R = v2, a small discrepancy 
from the proposed formula can be seen, which is very likely 
due to stronger finite-size effect around this point. 



R = \ 1 2 arccos A 

7T 



1 < A < 1. 



(60) 



When M ^ 0, R can be determined numerically by solv- 
ing the integjjal equations obtained from the Bethe-ansatz 
method.cltoca We set M at simple fractions 0, 1/5, 1/4, 
and 1/2 so that we can examine the dependence on the 
system size L. As in the critical dimer models studied in 
the previous section, the entropy S well obeys the scaling 
form S ~ aL + So + b/L. The subleading constant So 
obtained by fitting the data is plotted as a function of R 
in Fig. ||, which shows a remarkable agreement with 



So =\ogR- 



1 



(61) 



When increasing A at M = 0, the system undergoes a 
Kosterlitz-Thouless transition at A = 1 from the critical 
phase to a massive Neel phase with doubly-degenerate 
ground states. In a finite-size system, the double degen- 
eracy in the Neel phase is slightly split, and the ground 
state can be approximated by a macroscopic superposi- 
tion of ordered states. When A — > oo, such a state is 
given by 

lfl) = ^(IU...) + Ut...». (62) 

This st ate g ives S = So = log 2. As in the discussion 
of Sec. [IIP , one can expect that quantum fluctuations 
around the state ( |62| ) occurring in A < oo produce only 
extensive contributions and that the constant So = log 2 
is stable in the massive phase A > 1. Our numerical 



result for So is presented in Fig. £)[ The data show devi- 
ation from log 2 when decreasing A, but it is likely due 
to finite-size effects. The peak and dip seen in the fig- 
ure move to the left as we use larger L's for extracting 
So- In the thermodynamic limit, we expect a jump from 
So = logv2 — 1/2 to So = log 2 at the transition point 
A = 1. 

Finally, we note that the XXZ chain with A =_L/2 
and h = corresponds to the so-called ice modelJilEJ 
where all the configurations satisfying the ice rule (two- 
in and two-out around every vertex) occur with equal 
probabilities. The RK state built from the ice model has 
been studied for ths-spin and fermionic models on the 
checkerboard latticeB3'Eil The result in this section shows 
that the half-cylinder entanglement e ntro py of this state 
has a subleading constant So = log \/4/3 — 1/2. 



V. ISING CHAIN IN A TRANSVERSE FIELD 

As an example showing ac~ 1/2 critical point, in this 
section, we study an Ising chain in a transverse field: 



L-l 



L-l 



(63) 



3=0 



3=0 



This model is related to two types of 2D-classical models 
depending on which basis we work with.c^lln the a x basis, 
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XXZ chain : phase transition for the subleading term in the entropy 



Exact diag, fit L=20,22,...,32 
Exact diag, fit L=12,14,...,22 
Expected limit for S„ 
log 2 




FIG. 9: (color online) The constant part So of the entropy 
and the phase transition in the XXZ chain at M = 0. In the 
thermodynamical limit, we expect So to be So = logi?(A) — 
1/2 for — 1 < A < 1, and So = log2 for A > 1. 



we have a 2D Ising model: 



(64) 



where jf runs over all the nearest-neighbor pairs of 
sites on the square lattice. This model shows a low- 
temperature ordered phase and a high-temperature para- 
magnetic phase. On the other hand, in the a z basis, we 
have an eight-vertex model of special type. The spins a z 
are placed on the bonds of the square lattice and satisfy 
local constraints; the product of four spins around each 
site must be even: 



n 



= +1. 



(65) 



Then the four spins can take 8 states out of 2 4 possibili- 
ties, hence the naming, eight-vertex. The energy is given 
by 



(66) 



It is useful to introduce a loop representation of the con- 
figurations. We regard the lowest-energy state (a z = +1 
for all j) as the "vacuum," and place a loop element on 
every bond j with u| = —1. Then only closed loops 
are formed because of the local constraints (IgBT). Equa- 
tion ( |66|) means that the energy cost to generate loops is 
proportional to their total length. At low temperatures, 
the system contains only small loops and is dominated 
by the vacuum ("small-loop" phase). At high tempera- 
tures, the formations of large loops are allowed and the 
system gets disordered ("large-loop" phase). The corre- 
spondence among quantum and classical models is shown 
in Table |, together with our results for the entropy which 
we present below. 



Here we consider the Shannon entropies, and S^ x \ 
defined in the a z and a x bases respectively. These corre- 
spond to the half-cylinder entanglement entropies of the 
RK states built from the eight-vertex model and the 2D 
Ising model, respectively (notice that for the latter, one 
needs to modify the model slightl y in o rder to simplify 
the calculation, as presented in Sec. | IIB 2 j ) . The RK state 
constructed from the eight-vertex model ( p6|) in the large- 
loop phase is particularly interesting because it possesses 
topological order. Such a state has been studied-as the 
ground state of a quantum eight-vertex-4podeli3E-J (also 
known as an extended toric code modet-J). 

As is well known, the Hamiltonian ( p3| ) reduces to a 
fermionic quadratic form using the Jordan- Wigner trans- 
formation. It can then be diagonalized using the Bo- 
goliubov transformation (see Appendix |j). The weight 
|(j|<?)| 2 of each spin configuration \i) in the a z basis can 
be obtained by calculating a Pfaffian, and S^ Z \L) is com- 
puted numerically by summing over all the 2 L configura- 
tions. 

The scalings of S^(L) are shown in Fig. [n]. We again 
observe nice agreement with a linear scaling S^ Z \L) = 

aL + Sk both in the critical and massive cases. At the 



critical point (/i = 
the entropy to be 



1), we find the subleading constant in 



S { Q z) = -0.4387 ±0.0001. 



(67) 



Figure |TT| shows the constant part Sq as a function of 
the coupling constant p. Away from the critical point, 



O) 



in the 



the constant is stable at certain values: S t 

disordered phase (p < 1) and Sq = — log 2 in the fer- 
romagnetic phase (n > 1). Deviations from these val- 
ues near the critical point are likely due to finite-size ef- 
fects because they decay as we increase the system size. 
These values can be understood by considering two lim- 
its. In the limit p = 0, the ground-state wave function 
is \g) — |tt ■ ■ • T)zj an d the entropy is zero. In the 
limit [i — ► oo , the wave function is 



b) = -4=(ITT ... 1% + |U ... I}*) 



V2 



(68) 



V2 



= T ■ ■■ T 

C1 = T,1 0"L=T,1 



where all the configurations have an even number of up 
spins in the u z basis. Therefore, the entropy is 



sW(L) = (L-l)log2. 



(69) 



As explained in Sees. [IIID and [V, we expect that quan- 
tum fluctuations around these limits produce only ex- 
tensive contributions and keep the subleading constants 
stable. 

It is useful to interpret these results in terms of the 
eight-vertex RK state. When the temperature of the 
eight-vertex model is set to infinity, the corresponding 
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TABLE I: Correspondence between the Ising chain in a transverse field and related 2D classical models, and the results for 
the constant part of the entropy and S (x K 



Ising chain in a 
transverse field (163) 


disordered phase /i < 1 


c = 5 critical point fi = 1 


ordered phase (J, > 1 


Constant part 
of the entropy 


si z) = 

S ( X) = 


S ( Z) = -0.4387(1) 
= S ( z) + log 2 = 0.2544(1) 


S«=-log2 
=+log2 


Eight-vertex model (|r36"|) 


small-loop phase 
(low temperature) 




large-loop phase 
(high temperature) 


2D Ising model @ 


disordered phase 
(high temperature) 




ordered phase 
(low temperature) 



Ising chain in transverse field : entropy 



Subleading constant in the Ising chain in transverse field for general u. 







H=1 (critical) -| 




extrapolation: S =-0.4387 




u=0.5 (disordered) ■ ■■■)< 




extrapolation: S =-1e-07 






x /• x x 



FIG. 10: (color online) Entropy S^ z ' computed from the 
ground state of an Ising chain in a transverse field up to 
L — 36 spins. The data for /i = 1 (critical point) are well 
reproduced by 5 (z) ~ aL + S^ +8/L with S ( z) ~ -0.4387(1) 
(determined from a fit to the last three points L = 32, 34, 36). 
For /i = 0.5 in the disordered phase, the constant is very 
close to zero (a fit to the three points L — 18, 20, 22 gives 



?Wl 



< 10" 



RK state is an equal-amplitude superposition of all the 
loop configuration*. This is the ground state of Kitaev's 
toric code modelH Using the method of Ref. |l^, the 
half-cylinder entanglement entropy is shown to scale ex- 
actly as Eq. (]69|), and the constant part Sq rjr=n — log 2 
can be interpreted as the topological entropytftL3 associ- 
ated with Z2 topological order. Microscopically, — log 2 
comes from the fact that in any configuration, the loops 
cross an even number of time with the boundary separat- 
ing the two half-cylinders. The result in Fig. |ll| therefore 
demonstrates the stability of this topological entropy for 
the half-cylinder geometry in the entire large-loop phase. 
The jump in the entropy around fx = 1 is interpreted 
as a breakdown of topological order. We comment that 
the stability of topological entropy has also been stud- 
ied on the same wave function for the disk and annulus 



geometries in Refs. tta and 16 




FIG. 11: (color online) Subleading constant extracted 
from the entropy S^ Z '(L) in an Ising chain in a transverse 
field for different values of fi. 

Subleading constant in the entropy for the Ising chain as a function of p 

' S (uptoL=16) — | — 
S (up to L=26) -X- 
S (uptoL=36) 

0.4387(1) 

-log 2 




FIG. 12: (color online) Subleading constant Sq Z '(/3) extracted 
from the entropy S^ Z '(L;(3) in an Ising chain in a transverse 
field at the critical point fi = 1. The inverse temperature f3 
is introduced as explained in Sec. [ID. 



We move on to the entropy S^ x ' in the a x basis. It can 



13 



be related to £& z) 
transformation^: 



using the Kramers- Wannier duality 



Spin 1/2 XXZ chain : sub-leading constant in -log p 



(70) 



Here, <rj 
and crj 



-1 is 
Tak- 



by which H((j,) is related to Ti(l/n) 
identified with a domain wall between cum 
ing into account the two-to-one correspondence between 
o x and a z configurations, one can show 



S^Gu) = S«(l/M) + log2. 



(71) 



Hence we obtain the results summarized in Table |. Now 
we have a positive constant So = log 2 in the ordered 
phase fi > 1. This is a consequence of the macroscopic 
superposition of two ordered states, as discussed for the 
ordered phase of the XXZ chain in Sec. |v|. 

We have obtained two constants Sq Z = —0.4387(1) 
and S^ x) = S^ z) + log 2 = 0.2544(1) at the critical point, 
depending on the choice of basis. We expect that these 
might be generic constants characterizing the c = 1/2 
criticality, although at present we do not have any ana- 
lytical derivation of these numbers. 

We can also introduce a tem perat ure for the en- 
tropy as described in Sec. [II D| , The constant part 
Sq Z \(3) extracted by fitting S^ Z \L; 0) with a linear scal- 
ing changes rapidly around (3 = 2 (see Fig. FT3) . It seems 

( z) — 

reasonable to conjecture that Sq {(3) becomes a step 

( z) 

function in the thermodynamic limit : Sq {(3) = — log 2 

for p < 2 and S { z) (J3) = for p > 2. If confirmed, this re- 
sult would suggest that increasing P has a role similar to 
decreasing fx, i.e. taking the system away from its critical 
point. In the c = 1 case, So(P) was a smooth function 
of p (see Figs. [| fj] and ||). Therefore, it looks like intro- 
ducing P has a qualitatively different effect, depending 
on the nature of the critical theory. 

Note that there are several directions in which to ex- 
tend the Ising model. One possibility is to study the 
g-state Potts model or the RSOS models along the same 
lines. Another one is to view the Ising model as a spe- 
cial case (n = 1) of the dilute 0(n) loop model. In the 
loop model case, the are the probabilities that the 

equator of the sphere is run across by loops at positions 
/Xj. In that case, we would find a universal curve So(n) 
extending Eq. @. 



VI. SCALING OF THE LARGEST 
PROBABILITY 

In this section, we study the scaling of the largest prob- 
ability 



p := maxft = \ (io\g)\' 



(72) 



i.e., the weight of the most probable configuration i in 
the ID wave function \g). In terms of a 2D RK state, 
this corresponds to the largest eigenvalue of the reduced 



Exact diag M=0 + 
Exacl diag M=1/5 
Exact diag M=1/4 X 
Exact diag M=1/2 □ 
log[R] 




FIG. 13: (color online) The subleading constant 7 in the scal- 
ing of — logpo [see Eq. ([73])] extracted from the critical ground 
state of the XXZ chain in a magnetic field. The same setting 
as Fig. |I] is taken. 



density matrix pa of a half cylinder. Very similarly to the 
entropy S, we find that — logpo exhibits a linear scaling 
with L followed by a subleading universal constant: 

-logpo =aL + 7 + 0(1). (73) 

Below we evaluate 7 in some critical systems. 

A. c = 1 critical systems 

We first consider the XXZ chain in a magnetic field, 
Eq. d5S|), in the critical phase. We find that the largest 
probability po = |(«o|3)| 2 is attained by crystal states. 
For example, |» ) = |UU ■ ■ •> and ||tlt ■ ■ •) for M = 0, 
and |i ) = ITTT4TTT4- ■ ■), etc. for M = 1/2, independent 
of A (> —1). The constant 7 is extracted by fitting finite- 
size data with Eq. (73). As shown in Fig. |l3], we observe 
a simple relation: 



7 = logi?. 



(74) 



The same result can be shown exactly for the largest 
probability po{P) '■= maxpi(/3) in the Dyson-Gaudin gas 
when P = 2R 2 = 2A with A G N and A < see 
Appendix [d|. For general P, Eq . ( f 74|) can be numerically 
demonstrated as shown in Fig. |14|. 

In order to understand the connection between S and 
— logpo; it is useful to introduce the Renyi entropy: 



S (N) ._ 



^T los (? pf )' 



(75) 



where N is a real number. Then, S and — logpo corre- 
spond to the limits — > 1 and 00, respectively. Now 
we assume that the probability distribution is given by 
the Boltzmann weights Pi(P) = p { /Z{P) of the Dyson- 



Gaudin gas in the critical phase R 



The 
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Dyson-Gaudin gas : subleading constant in -log(p ) 
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FIG. 14: (color online) The subleading constant 7 in the scal- 
ing of — logpo(/9), where Po{f3) is the largest probability in 
the Dyson-Gaudin gas with density p = 1/3. The data are 
consistent with Eq. (|7J1. For some special values of 7? (see 
the text and App endix^). Gaudin's formula ( ^i|) can be used 
to show Eq. (|74|) exactly. 



In fact, Hsu et alB3 proposed the right hand side of 
Eq. ( [79]) as the expression of the entanglement entropy 
5 VN of a half cylinder. Their argument was based on a 
replica trick; the iV-th moment A4^ N ' := Tr = Y^iPf 
of the reduced density matrix was evaluated for inte- 
ger N > 2 and then an analytic continuation N — > 1 
was taken. According to their evaluation, the Renyi en- 
tropy S( N > = jPzi log M. ( N ' is expressed by the r.h.s. 
of Eq. ( |7S| ) for any integer N > 2, leading to an N- 
independent subleading constant logi?. On the other 
hand, we have obtained N -dependent constant (|77j). The 
two results for the subleading constant agree only in the 
limit N — > 00. We infer that this discrepancy comes 
from a difficulty in specifying boundary conditions in the 
argument of Hsu et alE3 In their argument, they took 
linear combinations of plural compactified fields, which 
could make the compactification conditions ambiguous. 
A more careful treatment of the compactification condi- 
tions and a derivation of Eq. ( |77| ) in boundary CFT are 
left as important open issues. 



Renyi entropy (|75[) can then be expressed as 



-1 
N - 1 



log 



Z(N0) 



(76) 



Recalling Eq. ( |55| ) for the non-extensive part C of Z((3), 
the subleading constant contribution to S' 1 ^ is given by 

f logi?- log N, VNR<R c = d, 

\ w=j logi?- jv3i log d, ^NR>R c = d. 

Both expressions give 7 = logi? in the limit TV — > 00. On 
the other hand, the former expression is consistent with 
So = logi? - 1/2 in the limit N -> 1. 

An alternative strategy to derive Eq. ( |74| ) is to adopt 
the 2D viewpoint of Sec. [fi], where the probability pq was 
related to classical partition functions: 



ZfZf 

logpo = - log l ° 10 



(78) 



Here. 



»0 



and Zf are partition functions on A and B 



with spins fixed in a state %q at their common boundary. 
Now we move on to the continuum limit described by the 
action (^). Recalling that io is given by a crystal state, 
the above boundary condition corresponds to locking the 
field at a certain cpastant at the boundary (Dirichlet 
boundary condition) Hence, we obtain 



logpo = - log 



■p-A 7 B 



(79) 



where D stands for Dipkhlet. This expression has been 
evaluated by Hsu et alBB and by Campos Venuti et alEB 
using boundary CFT. Their resul^Snfor the non-extensive 
part are consistent with Eq. (T74 



B. c = 1/2 critical system and beyond 



We next consider the Ising chain in a transverse field, 
Eq. j6^), at the c = 1/2 critical point fj, = 1. Figure [l5| 

shows the scaling of — logPo m the a z basis. Here, 
the largest probability is attained by the ferromagnetic 
configuration |io) = [ TTT - • -)z- In this case we have the 



following exact formula (see Sec. E3): 



(2) 
Po 



L/2 

n< 



(2j ~ 1)7T 

4L 



(80) 



An Euler-Maclaurin expansion shows that the subleading 
constant is = 0. The data in Figure [lj] is consistant 
with this result and gives 7^ < 10~ 6 . It can also be 
shown that remains zero away from the critical point. 
Moving to the a x basis using the same argument leading 
to Eq. @, we obtain 7^) = log 2. 

In a 2D viewpoint, po is related to a ratio of partition 
functions, as shown in Eq. (78). Since the boundary con- 
figurations, ah and bh, at the upper and lower edges of 
the cylinder were arbitrary, we can glue these edges by 
identifying ah and bh and integrating them out. Alge- 
braically, we consider 



Po 



(io\T 2h \i ) Z iQiQ (L,2h) 



Tr T 2h 



Z P (L,2h) 



(h 



(81) 



Here, the numerator Zi i is the partition function of a 
long cylinder with boundary configurations fixed in the 
same state io at both edges. The denominator Zp is the 
partition function of a torus. 

A similar quantity has been considered in a rational 
CFT context. Therein, the fixed boundary conditions 
imposed in the numerator of Eq. (81) are replaced by 
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Ising chain in transverse field - ferro. state probability 



[i=1 (critical) 
Fit : -log(p„)=0.110*L 



bility of the "+" state in CFT, hence 

= — logp + = log 2 



(86) 



in agreement with our calculation. 

These results can be extended to A rn RSOS models with 
central charge c < 1. The configuration p with highest 
probability is in this case 



Po 



m(m + 1) 



sin 



7T ^ 

sin I — 

1 / \m/ 



(87) 



VII. SUMMARY AND CONCLUSIONS 



FIG. 15: (color online) Scaling of — logpj) in the Ising model 
in a transverse field at the critical point, calculated in the a z 
basis. The subleading constant 7^' is very close to zero: a fit 
to the last three points L = 96, 98, 100 gives < 10" 6 . 



conformally invariant ones a. The associated probabil- 
ity p a can be evaluated as explained in Ref. |5J. In the 
notations of Ref. K4, the result is: 



Pa 



si 



(82) 



where tj) % a are certain structure constants characteristic of 
the model, and S\ is the identity matrix element of the 
S matrix implementing the modular transformation. In 
the simplest case of an A m SU(2)k model, the bound- 
ary fields correspond to the vertices of the A m Dynkin 
diagram. Let {d a } be the Perron-Frobenius eigenvector 
of the A m incidence matrix, normalized so that d\ = 1. 
Then, each d a is the so call ed quantum dimension of the 
state a, and d = y/Yla is called the total quantum 
dimension. In this case, ip a = d a /d and Si = l/d, hence 



(83) 



We study the case where the Dynkin diagram is A3 and 
the possible states are +, free, — with quantum dimen- 
sions 1, \/2, 1, respectively. The probabilities are, in the 
a x basis 



P-\ 



P- 



1 

2 ) Plree 



1. 



(84) 



A3 also describes the Ising model and this enables us 
to confirm the numerical results at the beginning of this 
subsection. The ferromagnetic state in the a z basis, | ttt 
. . . |)z, may be regarded as a paramagnetic state in the 
a x basis or the "free" state in CFT. Equation (|§^) then 
gives 



7 



0) 



logpfrcc = 0, 



(85) 



in agreement with our numerical result. The largest prob- 
ability po in the a x basis may be regarded as the proba- 



The starting point of this study was to introduce 
the Shannon entropy of a ID ground state wave func- 
tion, which measures quantum fluctuations occurring in 
a given basis. 

Like other entanglement measures, we have seen that 
the scaling behavior of this entropy is essentially con- 
trolled by the long-distance correlations. Using a trans- 
fer matrix approach, we showed that this entropy can 
also be interpreted as the entanglement entropy of a half 
cylinder for a suitably chosen 2D RK state. This corre- 
spondence allowed us to study the entanglement entropy 
of 2D wave functions using simpler ID systems, without 
the need to trace explicitly over the degrees of freedom 
sitting outside of the subsystem (a formidable task in 
2D). 

To unveil the generic scaling properties of the Shannon 
entropy of ID states (equal to the entanglement entropy 
of 2D states), we have studied several ID quantum sys- 
tems: (i) a discretized version of the Dyson gas/Calogero- 
Sutherland ground state wave function (relevant to 2D 
dimer RK states) in Sec. Ill, (ii) the spi n-i XXZ chain 
(relevant to six- vertex RK states) in Sec. [F^and (iii) the 
Ising chain in transverse field (relevant to 2D Ising RK 
states and 2D eight-vertex RK states) in Sec. [v| 

In both critical and massive systems, we found that 
this entropy is composed of an extensive part and a sub- 
leading constant So- There is no logarithmic contribution 
as anticipated before for half-cylinder geometry.E§Ej For 
Tomonaga-Luttinger liquids (cases (i) and (ii) above), 
described by a compactified boson with radius R, we 
showed numerically and analytically that So = log R — h 
(a result which differs from the recent prediction by Hsu 
et aLe3). Going back to the 2D entanglement entropy 
interpretation of this result, it implies that the usual 
RK states for dimers on the hexagonal or square lat- 
tice (with R = 1) have So = — 5. At present, we do 



not have a derivation for the value S, 



-0.4387 (or 

(x) 

S = —0.4387 + log 2 depending on the choice of the 
basis) found numerically for the Ising chain in transverse 
field at its c = 1/2 critical point. 

We introduced a temperature to extend the Shan- 



non and entanglement entropies in Sec. |II D| . This has 
different consequences depending on the nature of the 
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criticality. For TLLs, changing gives a natural way to 
tune the boson radius R while retaining the c = 1 critical- 
ity. This allowed us to "deform" a dimer RK problem (or 
a free fermion problem) and to derive the /independence 
of the entropy constant So in Sec. HL When (3 reaches 
a critical value, the system undergoes a phase transition 
to a crystal state, where the entropy constant takes a 
stable value So = logc? with d being the degeneracy of 
the ground states. Hence, this transition can be detected 
through S . In contrast, in the critical Ising chain with 
c = 1/2, the entropy constant So shows an abrupt change 
around (3 = 2. More generally, the /^-dependence might 
offer a valuable fingerprint for clarifying the nature of the 
undeformed case [3 = 2. 

We also considered the scaling properties of po, the 
probability of the most likely configuration in a critical 
ID state (Sec. VI). This quantity, which corresponds to 
the largest eigenvalue of the reduced density matrix in 
the 2D point of view, also contains a universal constant 
contribution 7 in critical states. We found numerically 
and analytically that 7 = log R for TLL states and 7^ = 
and 7 ( - :E ) = log 2 for the critical Ising chain. These are 
related to the probabilities associated with conformally 
invariant boundary conditions. 

The universal independence of So and 7 found in the 
present work can be used as a new simple way to deter- 
mine the boson radius R in a TLL through a ground- 
state structure. Similar universal i?-dependence was also 
found in the mutual information (double-interval entan- 
glement entropy) studied in Refs. [?| and ||. Notice that the 
mutual information is invariant under the transformation 
R — » 2/R while the present quantities are not. This dif- 
ference is related to the origins of the -R-dependence: it 
comes from certain boundary effects in the present case, 
while it comes from the special topologies of theJjiiemann 
surfaces in the case of the mutual informationM 

In fact, the transfer matrix approach gives access to 
all the eigenvalues {pi} of reduced density matrix pa of 
a half cylinder, (as a first application, the gap pi/po is 
computed in Appendix [5| for dimers on the hexagonal 
lattice). The present approach therefore provides a con- 
venient topi, to study the properties of the "entanglement 
spectrum'!^ in a 2D quantum state. 
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APPENDIX A: FREE FERMIONS AND WICK'S 
THEOREM 



The probabilities of Eq. ( p6| ) involve quantities such as 
(0|aiG2 • • • &2n|0), where the aj are linear combinations 
of fermion creation and annihilation operators. Wick's 



theorem then gives: 
(0|ai . . .a 2n |0) = (ai...a 2 „) 

%!<...< i n 

Vfc, ik < jk 
= Pf A, (Al) 

where e(a) is the signature of the permutation which 
transforms {1,2, . . . , 2n} into {h, 31,12, h, ■ ■ ■ in, 3n} ■ Pf 
denotes the Pfaffian. A is an antisymmetric 2n x 2n ma- 
trix given by : 



.4 





(aia 2 ) 



\aia 2 ) 




\aia 2r , 

(a2a2r, 



(A2) 



-(aia 2n ) -{a>2a2n) • ■■ 
The two following properties of Pfaffians are useful : 



Pf 



(Pf Af 

B 
B T 



detA 
±detB 



(A3) 
(A4) 



and allow fast numerical calculations using determinant 
routines. 



APPENDIX B: TRANSFER MATRIX FOR THE 
CLASSICAL DIMER MODEL ON THE 
HEXAGONAL LATTICE 

1. Transfer matrix as free fermions 

Here we consider a hexagonal lattice with periodic 
boundary conditions and an even number of columns L. 
The mapping onto free fermions is as follows (see Fig. fl6| ): 

We choose a convenient dimer configuration which we 
call a reference configuration. Any other dimer configu- 
ration (real configuration) will be compared to the refer- 
ence by superposition of the two. We define the particle 
locations as the vertical edges that are not occupied by 
a "real" dimer (only a reference one). Particles can jump 
from a vertical edge to another only if a real horizon- 
tal dimer connects the two. This mapping has several 
interesting properties : 

• The dimer configuration is totally determined by 
the trajectories of the particles. 

• Two particles cannot go to the same edge. There- 
fore, they obey a fermionic exclusion rule. This 
encodes the dimer hardcore constraint. 

• The number of fermions is conserved, so that the 
TM is block diagonal, each block corresponding to 
a fixed number of fermions. It should be remarked 
that this property would not hold on non-bipartite 
lattices (such as the triangular). 
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FIG. 16: (color online) Upper left : reference configuration. 
Upper right : real configuration. Below : Transition graph. 
Reference dimers are in blue. The fermions are living on the 
vertical edges of the lattice and are symbolized by red zigzag 
lines. The integers attached to each plaquette of the lattice 
form a height configuration associated to the dimer covering. 
When coarse-grained, these microscopic heights become the 
free fieldniwhich describe the long-distance properties of the 
system£jl£3 The heights can be constructed by fixing h = 
at some origin and then moving from plaquette to plaque- 
tte by turning clockwise 'around the sites of the even sublat- 
tice (marked with a black dot). The rule is the following: 
the height picks a contribution equal to +2 when crossing a 
dimer, and -1 otherwise. Since there is exactly one dimer 
touching each site, the height difference between two points 
does not depends on the chosen path on a simply connected 
domain. With periodic boundary conditions, the height is 
not single-valued. For example, when winding horizontally 
around the system, the height picks a contribution W x (also 
called winding number) equal to twice the number of verti- 
cal dimers crossed, minus the number of empty bonds. In- 
serting a fermionic world line going upward shifts the height 
by —3 by going from the left to the right, and thus changes 
W x — ► W x — 3. It is simple to check that the configurations 
with a fermion density equal to | have W x = and an average 
"slope" equal to zero. 



Fermionic representation and periodic boundary 
conditions 



Therefore, to keep the order of Eq. (Bl) : 

L-l 



{-l) n -\ 



u 

o 



3=0 



(B3) 



In the following, we will also need the L fermion operators 
in Fourier space : 



,.t 



1 L ~ 1 

rr Z^ e V 



They satisfy {ck ,c' k ,} = 6kk' , provided 
The set of wave- vectors is given by 

Z = 0,...,L-1 



ikL 



(B4) 



7r 2nl 

+ z + - 

ke { - n + — 



Z = 0,...,L-1 



n odd. 



3. Diagonalization of the Transfer matrix 

Each fermion can go to the left or to the right with 
equal amplitude. We number vertical edges in such a 
way that a fermion located on j can go to j or j + 1. T 
satisfies 



T|0) = |0) 
TclT- 1 = cl + c 



t 



So that 



TclT- 1 = \{k)c{ , A(fc) = 1 + e 



ik 



(B5) 
(B6) 

(B7) 



c k \0) is then eigenvector of T with eigenvalue X(k). In a 
similar manner : 

Tclcl . . . clj- 1 = A(fcO . . . A(fc„)c t fei . . . 4 n . (B8) 

Provided all the wave vectors are different, . . . c' k |0) 
is an eigenvector with eigenvalue A(fci) . . . A(fc„). The 
transfer matrix can also be expressed explicitly : 



T=l[(l + e ik c{c k 



(B9) 



A state of a row is determined by the number n of 
fermions and their positions < ai <,...< a n < L — 1. 
We can choose to represent such a state using second- 
quantized fermions creation operators : 

|i) = |ai J ...a n )=4 1 ...ct n |0). (Bl) 

As we want to use the translational invariance, we have 
to set 

|ai, . . . ,a n -x,L) = |0,ai, . . . ,a n _i). (B2) 



4. Largest eigenvalue and dominant eigenvector 

Since |A(fe)| > 1 for every k e [-2tt/3, 2tt/3] (see 
Fig. the eigenvalue with largest modulus in a given 
sector with n fermions is obtained by a product over the n 
nearest to wave vectors. Let us denote this eigenvalue 
by A„. Then, A max = Max {|A„| , 1 <n< L}. The 
eigenvalue with largest modulus is real and has approxi- 
mately all allowed k lying in the interval [— 2tt/3, 27r/3]. 
The dominant sector has therefore n ~ 2L/3 fermions. 
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FIG. 17: (color online) One-particle eigenvalue as a function 
of k 



It is easy to understand the fact that the "dominant" 
fermion density is 2/3 : it corresponds to flat height con- 
figurations (see Figjlq). 

If we denote by ft the set of wave-vectors that gives the 
largest eigenvalue then the dominant eigenvector is : 



Is) 



n* 



.t 



|0>. 



(BIO) 



Let us show what is Q in the simple case where L = 6p. 
We have to distinguish between the even and odd sectors 




FIG. 18: (color online) Upper left : reference configuration. 
Upper right : chosen configuration. Below : the fermions are 
living on the vertical edges of the lattice and are symbolized 
by red zigzag lines. Edges are numbered from to L — 1. The 
integers attached to each plaquette form a height configura- 
tion associated to the dimer covering. The rule is very similar 
to that of the honeycomb lattice : turning clockwise around 
the sites of the even sublattice (marked with black dots), the 
height picks a contribution equal to +3 when crossing a dimer, 
— 1 otherwise, h = is fixed at some origin. For a more de- 
tailed presentation, see Ref. 



5p-l 

= Max{A 2n ,} = Y[ A ( -tt + tt 

i=p 

Alt = Max{A 2 „, +1 } = J] A ( -tt + ^ 

i= P 



21 + 1 
6p 

nl 

3p 



(Bll) 



Here, Amlx > Amax because Euler-Maclaurin expansion 
gives log A^L- log A^L = t$ + o(1/p). Therefore the 
leading eigenvalue corresponds to Ap = |i fermions and 

21 + 1 



6p 



P < I < 5p - 1 



(B12) 



We get a Vandermonde determinant and pi simplifies into 
Eq. ©. 



APPENDIX C: TRANSFER MATRIX FOR THE 
CLASSICAL DIMER MODEL ON THE SQUARE 
LATTICE 



5. Probability of a given configuration 



Free fermions 



The dominant eigenvector has n = 2L/3 fermions. A 
configuration i is represented by 

N> = 4... c tjo>, 



(B13) 



and will have a probability (we use Eq. (Al) and 
Eq. @) : 

2 



Pi 



"ai c \ 1 



T ) det(e-^V) ij ,| 



(B14) 
(B15) 



We consider a square lattice with periodic boundary 
conditions and an even number of columns L. The map- 
ping is similar to that of the hexagonal lattice. The ref- 
erence configuration is shown in Fig. Here, a fermion 
will be defined as an even vertical edge occupied only by 
a reference dimer, or an odd vertical edge occupied only 
by a real dimer. It can go to the left, straight ahead, or 
to the right. We introduce a shift in the numbering, so 
that a fermion located on site j can go to j, j + 1, or 
3 + 2. 



19 



2. Diagonalization of the transfer matrix 

As for the honeycomb case, T is block-diagonal and 
invariant by translation. It also satisfies 



Tc^T 

1 ^9,J_l 1 



T\0) = |0) 

t T-l 



C 2j+2 + C 2j + 1 + C 2j 
J 



"2J + 1- 1 - L '2j+2- 

As usual we also define Fourier space fermions 



4 

- OA- 
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2; 



-ifc(2j+l)„t 



-2j + l' 



(CI) 
(C2) 
(C3) 



(C4) 
(C5) 



with 



2 

kel - 



7T 


2?r/ 




" ~L 




2ttI 


2 " 





1 = 0,, 
I = 0, 



- 1 
- 1 



n even 



n odd. 



The transfer matrix acts on them in the following way 



1 c 0k 2 



(1 + e 2lk )c! 



Ok 



e c 



Ik 



t n 1 



7V T 

1 c lk 2 



ikj 
e c 0k- 



(C6) 
(C7) 



To diagonalize T, it is therefore sufficient to diagonalize 
a 2 x 2 matrix : 




|A + (fc)|- 
|A_(fc)|- 



7T/4 



7T/2 



FIG. 19: (color online) One-particle eigenvalues A + and A_ 
as a function of k 



3. Largest eigenvalue and dominant eigenvector 



We assume for s imp licity that L is a multiple of 4. 
Noticing (see Fig. that Vfc, |A+(fc)| > 1 and also 
|A-(fc)| < 1, we can deduce that only the A+( fc) will con- 
tribute to the largest eigenvalue. In Ref. it is shown 
that 



I A+(fc), 



(C15) 



ken 



where = {-§ + f + fl | I = 0...L/2 - 1}. The 
leading sector has an even number of fermions L/2 and 
the dominant eigenvector will be given by 



M 



1 + e 2lk e ik 



(C8) 



\g) 



n & u i°>- 

\ken / 



(C16) 



If one sets 



then 



tan 6k 
A±(ft) 



\/l + cos 2 fc — cos fc 
e ifc (cos k ± \/i + cos 2 k 
cos 9 k cl k + smQ k c\ k 



= - sin 9kc' ok + cos 9kc\ 



k> 



Tb^T- 1 = X+(k)b\ 



J +k 



+k 



X-(k)bl k , 



(C9) 
(CIO) 

(Cll) 
(C12) 



(C13) 
(C14) 



which gives us all the eigenvalues and eigenvectors of T. 
It is also possible to express explicitly T : 



4. Probability of a given configuration 

The dominant eigenvector has n = L/2 fermions. A 
configuration i is represented by: 

W=4---4J0>, (C17) 

and will have a probability (using Eq. (Al) and also 
Eq. @): 

2 



Pi 



...c ai b 

|det m^- 



a i u +,k 
2 



10) 



where 



cos6» fe , e i£ W 
sin Qk.,e laik i' 



otj even 
ay odd 



(C18) 
(C19) 

(C20) 



This determinant is slightly more complicated than on 
the honeycomb lattice, and cannot be further simplified. 
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APPENDIX D: 2D COULOMB GAS ON A 
CIRCLE 

We consider a Gaudin model with n charges dispatched 
on a circle with L sites. 



1. Scaling of the ground state 

We study the special case where L/n G N and we set 
p = n/L. We denote by po(P) the probability associ- 
ated with the ground-state configuration. It corresponds 
to the case where the distance between each charges is 
maximal. Therefore it is obtained when all charges lie on 
the edges of a polygon. Hence, 

P0(P) = — | L- nf3 ' 2 \{ | e 2^/«„ e 2,W«|/3. 

Zn \P) l<k<l<n 

Using the formula YVi=i (l — e 2l W«) = n , we get : 



1/3/2 



zi N \p) 



(Dl) 



In the special case where (3 is an even integer and (3 < 
it is possible to use Eq. ( j33[) and the subleading 
term of — logpo(00 can easily be found : 



1 (3 

logp (/3) = an + - log-. 



(D2) 



This result for the universal part of the probability, 7 = 
I log ^ = logi?, does not depend on L nor n (no finite- 
size effects). 



Q n {x) and P n {x) are of the same degree, and share the 
same zeros : they have to be proportional. P„(0) = 1 
and Q n (0) = n yields P„(x) = ±Q n (x). Using Eq. (|d|) 
and Moivre's formula, we get : 



[(n-l)/2] 



P n {x)= l - Y. M)^ 1 ; 



2/v 



fc=0 



Therefore P n (— ) reduces in the limit n — > 00 to: 



P n { — ) = — sin(Trp). 

71 7T/9 



Finally, 



Pi(g) 



— sin(7rp) 
np 



(D6) 



(D7) 



(D8) 



So that there is a finite gap in the thermodynamical limit: 

" 1 



AE = E\ — E 



log 



Tip 



■ sm(irp) 



(D9) 



This calculation can easily be extended (in the thermo- 
dynamical limit) to any configuration deduced from the 
ground state by moving a finite number of particle. For 
the corresponding RK wave-function, AE gives an infor- 
mation about the first gap of the reduced density matrix 
(entanglement spectrum) . 



APPENDIX E: GROUND STATE OF THE ISING 
CHAIN IN A TRANSVERSE FIELD 



1. Diagonalization 



2. Existence of a gap in the thermodynamical limit 

The first excitation is obtained from the ground-state 
configuration by moving one particle to the next site, 
while keeping all the others in place. The associated 
probability pi (j3) will be given by : 



Pi 03) 
Po(P) 

In the limit n — > 

Pi 03) 
Po09) 



n 

1=1 



pj/n 



(D3) 



-00 it is possible to expand the sin : 



no 

1=1 



cot(?7r/n) 







(D4) 



We then consider P n (x) = Yl'iLi [1 — x cot(lir /ri)), The 
trick to calculate P n (x) is to introduce another polyno- 
mial Q n {x) which satisfies : 



Q n (tani) 



sin nt 



(sin t) cos" 



(D5) 



We consider the Hamiltonian of an Ising chain in a 
transverse field with an even number of sites L. 



L-l L-l 

Using a Jordan- Wigner transformation 

ju / j-i 



(El) 



_ 05 + 10 



3 - = c]exp [inJ24 c i] ( E2 ) 



1=0 



0] = 2c)c 3 1, 



(E3) 



TL is rewritten as 

L-l 



L-2 



H = - ^(2c] Cj --l)- M ^(cU Cj )(ct +1 +c 3 - +1 ) 
3=0 3=0 
+ //( C Li-c i _ 1 )(4 + co)e" Ar , (E4) 

where j\f = ^2^=0 c]cj is the fermion number operator. 
Since V = II^o 1 c| = ±1 is a conserved quantity, H 
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may be separately diagonalized in two sectors. Here we 
are only interested in the ground-state of the chain. In 
the basis of the eigenstates of the cr|, all off-diagonal 
elements are negative, and it lies in the sector V = +1 
(Perron- Frobenius theorem). Using V = exp(iirj\f) and 
the last term of Eq. (FA), we see that in this sector, 
one has to keep configurations with an even number of 
fermions satisfying anti-periodic boundary conditions 

4 = -4- (E5) 

To take advantage of the translational invariance, we in- 
troduce Fourier-space fermions 



L-l 



c f - 



—iki t 

e J c] 



(EC) 



3=0 



where k € {{21 + 1)£ | -L/2<1< L/2 - 1} are the L 
wave-vectors. The Hamiltonian becomes 

H= L - ^ 2(1 + yucos k)c\ck 

k 

+ M X] (i sin kc\cl k - i sin kc- k c k ) . (E7) 



This expression can be diagonalized using a Bogoliubov 
transformation 



-k> 



c\ = cos0- h - »-««-i-t 
0-k = —6ki 

provided the following condition is satisfied 

usinfc 

tan2(9 fc = — - . 

1 + /i cos k 

We then obtain : 

k 

e(k) 



(E8) 
(E9) 



(E10) 



(Ell) 



(1 + ficosk) cos26k + /isin/csin20fc.(E12) 

z, 

In the following, we want the vacuum of the bk to be the 
ground-state of H, which is true only if e{k) > Vfc. One 
also has to take into account the indetermination of 9 
modulo 7r. There are two cases : 

• // < 1 : Here 1 + /i cos k is always positive. We 
choose 



1 



arctan 



/' 



sin fc 



2 \1 + ji cos k > 

and the energy spectrum is given by : 



e(k) = 2^1 + 2^ cos fc 



(E13) 



(E14) 



• fj, > 1 : Here, one has to be careful because 
1 + /i cos k can vanis h and change sign. A generic 
solution of Eq. (E10) is 



1 / fi sin k \ 7r 

6» fe = - arctan — + -q k , g-fc G Z. (E15) 

2 \ 1 + [i cos k / 2 

The eigenenergies are given by : 



e(k) = (-l)*sgn(l + /iC0sfc)2Vl + 2^icosfc + p? . 

(E16) 

1 + yucosfc changes sign at k = ±k c = 
± arccos(— l/n). A possible choice is therefore : 

, -k c < k < k c ■ (E17) 

1 , k ^* k^ 

2. Probability of a given configuration 

Since |0) is the ground state of the chain, the proba- 
bility of each configuration i is (in the u z basis) 



(E18) 



where Pj (resp. Pj) is the projector onto the | |)| (resp. 
| |)|) state : 



PJ 



(E19) 



Using Wick's theorem, pi reduces to a Pfaffian. To com- 
pute it, we need to calculate four types of contractions : 
(cjCj>), (cjCj,), (CjCj,) and {cjCj'). This can be done by 
expressing back the Jordan- Wigner fermions in terms of 
the Bogoliubov fermions. 



-^cos 2 fc cos[fc(j -j')} (E20) 

k 

I^sin^cos^J-j')] (E21) 

k 

- ^2 sin k cos 6 k sin [k(f - j)} (E22) 

k 

- ^ sin B k cos 9 k sin [k{j - j')} . (E23) 



T /I 

If we write a generic projector as = a 2 j-ia 2 j, where 
a is either core', then 



(c t c t . 



Pi = Pf i 



<i,j<2L 



(E24) 



Notice that it is also possible to compute when L is 
odd. The only difference is that the fermion number op- 
erator satisfies 



cxp(injV) = —V. 



(E25) 
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Therefore, one has to take periodic boundary conditions 
c l = c 0' an d on ly keep configurations with an odd num- 
ber of fermions. The wave-vectors are now in the set 
{2ln/L | - L/2 < I < L/2 - 1}. Since the dispersion 
relation e(k) is minimum for k = —n, the ground-state 
wave function is 

|6)=6ljO), (E26) 
and the probabilities 

Pl = (0\Pl /l P} n . . -Pl /l \0) (E27) 
will also be given by Pfaffians. 

3. Configuration with highest probability 

We study the case where L is even. The configuration 
with highest probability is attained by the ferromagnetic 
configuration 

Ko> = | Tit •••}*• (E28) 

Defining 

L-l 

P=Y[c] Cj , (E29) 



Po is expressed as 



Po = <0|P|0). (E30) 



P is a projector onto a state with L fermions. 



Q = l[c{c k (E31) 

fc 



is also a projector onto a state with L fermions. Since 
there is only one state of the Hilbert space with L 
fermions, P and Q are in fact identical. Using Eq. (|Egj), 
we therefore get : 



Po = Y[cos6 k , (E32) 

fc 



where 6 k is given by Eq. (||l|) or Eq. (||l|). At the 
critical point {fi = 1), 9k = fc/4. 
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